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Abstract 

A  number  of  techniques  have  been  developed  in  recent  years  for  the  analysis  and  design 
of  controllers  which  are  robust  with  respect  to  structured  complex  uncertainty.  In  particular 
the  complex  p  synthesis  procedure  has  been  successfully  applied  to  a  number  of  engineering 
problems.  However  the  presence  of  real  parametric  uncertainty  in  the  problem  description 
substantially  complicates  matters,  so  that  standard  complex  n  synthesis  techniques  are  no  longer 
adequate.  In  this  paper  we  develop  a  procedure  to  tackle  the  mixed  (real  and  complex)  /i 
synthesis  problem.  This  procedure  involves  a  “D,G-K  iteration”  between  computing  the  mixed 
H  upper  bound  and  solving  an  'Hca  optimal  control  problem,  and  has  guaranteed  convergence  to 
a  local  minimum  of  the  (nonconvex)  problem.  The  procedure  has  been  implemented  in  software, 
and  several  controller  designs  are  compared  with  the  corresponding  complex  p  synthesis  designs. 


1  Introduction 

Many  control  engineering  problems  motivate  the  consideration  of  robustness  problems  involving 
both  real  and  complex  uncertainties.  These  mixed  uncertainty  descriptions  offer  the  control  engineer 
the  possibility  of  more  accurately  modeling  given  physical  systems,  and  hence  using  less  conservative 
analysis  tools,  to  produce  control  system  designs  with  enhanced  performance.  The  price  to  be  paid 
for  this  added  sophistication  to  the  analysis  and  design  procedure  is  an  increase  in  the  complexity 
of  the  associated  mathematical  problems  to  be  solved.  A  number  of  practical  analysis  tools  have 
been  developed  for  mixed  problems,  and  computational  software  is  now  available  for  the  general 
mixed  fi  analysis  problem  [YND92,  BDG+91].  In  this  paper  we  will  be  concerned  with  the  general 
mixed  n  synthesis  problem. 

The  problem  of  synthesising  a  controller  which  is  (optimally)  robust  to  structured  mixed  un¬ 
certainty  is  very  difficult,  since  the  associated  optimization  problem  is  not  convex.  Furthermore 
it  seems  intuitively  clear  that  the  synthesis  problem  is  at  least  as  hard  as  the  analysis  problem, 
which  is  known  to  be  NP  hard  [RP,  BYDM].  Some  exact  solutions  have  been  presented  for  special 
cases  of  the  synthesis  problem  (see  [RM92]  for  example,  which  reduces  the  rank  one  \i  synthesis 
problem  to  a  convex  optimization  problem),  but  these  are  all  cases  for  which  the  analysis  problem 
also  simplifies  considerably.  As  yet  there  is  no  globally  optimal  solution  to  the  general  synthesis 
problem  (even  in  the  purely  complex  case),  and  no  indication  that  one  will  be  forthcoming  in  the 
foreseeable  future. 

Nevertheless  the  (complex)  ^-synthesis  procedure  first  outlined  in  [Doy85]  has  been  successfully 
applied  to  a  large  number  of  engineering  problems  (see  [BYD92]  for  example).  This  procedure 
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involves  a  “ D-K  iteration”  between  computing  the  p  upper  bound,  and  solving  for  an  Hrj0  (sub) 
optimal  controller  (both  of  which  are  convex  problems).  This  procedure,  which  was  developed 
for  p  problems  involving  only  complex  blocks,  does  not  guarantee  to  find  the  globally  ^-optimal 
controller  (since  the  problem  is  not  jointly  convex  in  D  and  K ),  but  has  often  been  found  to  work 
well  in  practice.  The  approach  taken  here  to  the  mixed  p  synthesis  problem  is  to  extend  the  above 
procedure  to  the  mixed  case,  by  exploiting  some  new  analysis  tools  recently  developed  for  the  mixed 
p  upper  bound.  In  this  way  a  “ D,G-K  iteration”  procedure  is  developed,  which  finds  a  controller 
that  stabilizes  the  nominal  system,  and  attempts  to  minimize  the  peak  value  (across  frequency)  of 
mixed  p. 

Prior  to  tackling  the  mixed  p  synthesis  problem,  we  need  to  develop  some  State  Space  factor¬ 
ization  machinery,  and  this  is  the  subject  of  section  3.  A  number  of  factorizations  are  presented, 
which  are  required  in  the  mixed  p  synthesis  procedure  to  essentially  produce  stable  system  realiza¬ 
tions  of  the  scaling  matrices.  State  Space  formulae  are  presented  for  all  the  required  factorizations. 
This  machinery  is  then  used  in  section  4  to  develop  the  “D,G-K  iteration”  procedure  for  mixed 
p  synthesis.  The  procedure  utilizes  a  combination  of  Hco  optimal  control,  mixed  p  upper  bound 
computation,  and  the  above  factorization  machinery  to  yield  an  algorithm  for  mixed  p  synthesis. 
The  algorithm  has  guaranteed  convergence  to  a  local  minimum  of  the  problem  (as  in  the  complex 
case).  Full  details  of  the  procedure,  including  appropriate  routines  for  fitting  the  scaling  matrices, 
are  given  in  section  4.  This  procedure  has  been  implemented  in  software,  and  in  section  5  the 
procedure  is  applied  to  several  example  problems  involving  real  parametric  uncertainty.  It  is  seen 
that  the  new  procedure  is  able  to  exploit  the  phase  information  in  the  real  uncertainties,  to  produce 
a  controller  with  a  lower  value  for  mixed  p  (and  hence  better  robust  performance)  than  one  could 
obtain  with  standard  complex  p  synthesis  techniques. 

2  Notation  and  Definitions 

The  notation  used  here  is  fairly  standard,  and  is  taken  from  [You93].  For  any  square  complex  matrix 
M  we  denote  the  transpose  by  MT,  and  the  complex  conjugate  transpose  by  M* .  The  largest 
singular  value  and  the  structured  singular  value  are  denoted  by  W(M)  and  p,>c ( M )  respectively. 
For  a  Hermitian  matrix  M,  then  A (M)  denotes  the  largest  (real)  eigenvalue.  We  denote  the  k  X  k 
identity  matrix  and  zero  matrix  by  Ik  and  0k  respectively.  Occasionally  we  will  drop  the  subscripts 
from  I  and  0,  whence  they  denote  identity  and  zero  matrices  respectively,  of  the  appropriate  size. 

We  define  TZM  to  be  the  space  of  real-rational  proper  transfer  matrices.  For  a  transfer  matrix 
P  £  TZM,  then  we  denote  its  State  Space  representation  by 


and  we  define  P~(s)  =  PT(—s).  We  define  TZCm  to  be  the  subspace  of  TZM  with  elements  analytic 
on  Re(s )  =  0,  and  TZTioo  to  be  the  subspace  of  TZM  with  elements  analytic  in  Re(s )  >  0.  Given  a 
transfer  matrix  P  6  TZH^ ,  then  we  denote  its  infinity  norm  by 

||P||oo  =  sup  cr(P(jw))  (2) 

The  definition  of  p  is  dependent  upon  the  underlying  block  structure  of  the  uncertainties,  which 
is  defined  as  follows.  Suppose  we  have  a  matrix  M  £  <C"Xn  and  three  non-negative  integers  mr,  mc, 
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and  mc  (with  m  :=  mr  +  mc  +  mc  <  n )  which  specify  the  number  of  uncertainty  blocks  of  each 
type.  Then  the  block  structure  K,(mr,mc,mc )  is  an  m-tuple  of  positive  integers 


kmr ,  kmr+i  ,  .  .  .  ,  &mr+mc  5  ^mr+rac+l  ?  •  •  •  j  km) 


(3) 


This  m-tuple  specifies  the  dimensions  of  the  perturbation  blocks,  and  we  require  Ya=i  k{  =  n 
in  order  that  these  dimensions  are  compatible  with  M.  This  determines  the  set  of  allowable 
perturbations,  namely  define 


Xjc  =  { A  =  block  diag  (6[Ikl 6rmrIkmr ,  6[Ik 


ntf  *fwc  ' 


aC  aC  n 

,  ,  .  .  . ,  L\m  ) 


6C 

mr  +  l  >  '  "  •  > 

S-  G  R,S‘  £  C,  Af  G  Cfcm’-+mc+iX*’nr+mo+i}  (4) 


Note  that  XK  C  Cnxn  and  that  this  block  structure  is  sufficiently  general  to  allow  for  (any  com¬ 
bination  of)  repeated  real  scalars,  repeated  complex  scalars,  and  full  complex  blocks.  The  purely 
complex  case  corresponds  to  mr  =  0,  and  the  purely  real  case  to  mc  =  mc  =  0. 

Note  also  that  all  the  results  which  follow  are  easily  generalized  to  the  case  where  the  full  complex 
blocks  need  not  be  square,  and  the  blocks  may  come  in  any  order.  We  make  these  restrictions  in 
(4)  purely  for  notational  convenience. 


Definition  1  ([Doy82])  The  structured  singular  value,  of  a  matrix  M  G  <C"xn  with  re¬ 

spect  to  a  block  structure  K,(mr,mc,mc)  is  defined  as 


=  ftfiip  (cr(^)  :  det  (I  —  AM)  =  0} 


(5) 


with  nic(M)  =  0  if  no  A  £  X>c  solves  det(/  —  AM)  =  0. 

In  order  to  develop  the  relevant  theory  we  define  the  following  sets  of  block  diagonal  scaling  matrices 
(which,  like  p  itself,  are  dependent  on  the  underlying  block  structure). 


VK-{  block  diag  (D i, . . . ,  Dmr+mc ,  ^Amr  +  mc  + 1  >  •  •  •  >  dmCIkm  )  : 

Di  G  CkiXhi , det(A)  ^  0, dt  G  C, d{  ±  0}  (6) 

Qk  =  {block  diag (Gu...,Gmr,Okmr+1,...,Okm)  :  G,  =  G*  G  Cfc*xfc*}  (7) 


3  State  Space  Factorization  Theory 

Prior  to  considering  the  mixed  p  synthesis  problem  we  need  to  develop  some  State  Space  factor¬ 
ization  machinery.  We  will  develop  State  Space  formulae  for  specific  coprime,  inner-outer,  and 
spectral  factorizations  that  will  be  needed  for  the  mixed  p  synthesis  problem.  These  factorizations 
are  extensions  of  standard  results  (see  [ZDG]  for  example),  and  as  such  are  of  some  interest  in  their 
own  right. 

3.1  Basic  Concepts 

We  begin  with  a  brief  review  of  some  basic  concepts  about  transfer  matrices.  For  a  more  detailed 
exposition  on  this  subject  see  [ZDG,  Fra87]. 
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Definition  2  Suppose  we  have  a  transfer  matrix  D  £  71  Ad.  Then  D  is  said  to  be  inner  if  D  £ 
TlHoo  and  D~D  =  I,  and  co-inner  if  D  £  TTHoo  and  DD~  =  I. 

If  we  are  interested  in  square  transfer  matrices,  then  these  notions  can  be  simplified  somewhat. 

Lemma  1  Suppose  we  have  a  square  transfer  matrix  D  £  1ZM .  Then  the  following  conditions  are 
equivalent: 

1.  D~D  =  I. 

2.  DD~  =  I. 

3.  D  is  invertible  and  D-1~D_1  =  I. 

4-  D  is  invertible  and  D-1D_1~  =  I. 

Proof:  Follows  from  straightforward  manipulations.  □ 

Lemma  2  Suppose  we  have  a  square  transfer  matrix  D  £  IZAd .  Then  D  is  inner  if  and  only  if  D 
is  co-inner. 

Proof:  Apply  lemma  1  to  D  £  TZTi^,.  □ 

A  closely  related  concept  is  that  of  an  all  pass  transfer  matrix. 

Definition  3  Given  a  square  transfer  matrix  D  £  71  Ad,  then  D  is  said  to  be  all  pass  if  D  £  TZCco 
and  D~D  =  /. 

Note  from  lemma  1  that  we  could  equivalently  require  D  £  IZCoo  and  DD~  =  I.  For  square 
matrices  then  an  inner  (or  co-inner)  transfer  matrix  is  exactly  a  stable  all  pass  transfer  matrix. 

Lemma  3  Suppose  we  have  a  square  transfer  matrix  D  £  TIM.  Then  D  is  all  pass  if  and  only  if 
D  is  invertible  and  D  1  is  all  pass. 

Proof:  ( — ►)  Suppose  D  is  all  pass,  and  hence  D  £  TIC From  lemma  1  we  immediately  find 
that  D  is  invertible  and  D-1~D-1  =  I.  But  now  this  implies  D-1  =  D~,  and  hence  D-1  £  TIC 00, 
which  together  with  the  earlier  result  implies  D-1  is  all  pass. 

(< — )  Suppose  D  is  invertible  and  D"1  is  all  pass.  Applying  the  above  proof  we  obtain  D  =  (D-1)-1 
is  all  pass.  □ 

The  term  all  pass  is  motivated  by  the  fact  that  at  any  given  frequency  the  gain  matrix  of  an  all 
pass  transfer  matrix  is  unitary,  or,  in  other  words,  has  unity  magnitude  gain  in  all  directions.  To 
be  more  explicit,  we  have  the  following  well  known  results  (stated  without  proof). 

Lemma  4  Suppose  Ni,N2  are  all  pass  transfer  matrices,  and  D  £  TZM.  Then  for  any  w  £  R  we 
have 

cF(A1(jw)D(jw)A2(ju;))  =  a(D(jw)) 

Lemma  5  Suppose  Ni,N2  are  square  inner  transfer  matrices,  and  D  £  7THrM .  Then  we  have 
iViDiV2  £  TZTtoa  and 

HAjDiV.IU  =  ||D  Hoo 
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Thus  the  maximum  singular  value  frequency  response  of  a  system  is  invariant  to  an  all  pass  transfer 
matrix  (and  hence  the  Hoo  norm  is  invariant  to  a  square  inner  transfer  matrix).  This  invariance 
is  the  reason  we  are  interested  in  these  concepts.  The  p  upper  bound  is  based  on  computing  the 
maximum  singular  value  of  a  scaled  matrix,  and  it  will  be  seen  later  that  the  synthesis  problem  can 
be  tackled  by  employing  a  scaled  Hoo  norm  optimization  procedure,  which  will  rely  on  the  above 
invariance  properties. 

The  remainder  of  this  paper  will  make  use  of  the  following  notions.  The  term  stable  refers  to 
a  system  (matrix)  with  all  its  poles  (eigenvalues)  in  the  open  left  half-plane,  and  marginally  stable 
refers  to  a  system  (matrix)  with  all  its  poles  (eigenvalues)  in  the  closed  left  half-plane.  The  term 
antistable  refers  to  a  system  (matrix)  with  all  its  poles  (eigenvalues)  in  the  open  right  half-plane. 

Additionally  we  will  make  use  of  some  fairly  standard  notation  and  results  regarding  the  solution 
and  properties  of  Riccati  equations.  We  refer  the  reader  to  [DGKF89]  for  a  concise  review  of  this 
area,  and  [ZDG]  for  a  more  detailed  exposition. 


3.2  Coprime  Factorization 

Now  we  wish  to  develop  certain  transfer  matrix  factorizations,  the  first  of  which  is  a  coprime 
factorization.  Essentially  this  amounts  to  splitting  a  transfer  matrix  into  stable  numerator  and 
denominator  matrices,  without  allowing  stable  common  factors  between  them.  We  will  use  this 
concept  to  factor  out  the  antistable  part  of  the  original  transfer  matrix. 


Definition  4  Suppose  we  have  two  transfer  matrices  N,M  £  TZM.  Then  N,M  are  said  to  be 
right  coprime  (over  IZTioo)  if  N,M  £  TTH^,  they  have  the  same  number  of  columns ,  and 
left  invertible  in  TTHrj0 . 

Given  a  transfer  matrix  D  €  TZM,  then  a  right  coprime  factorization  of  D  (if  it  exists)  consists 
of  two  right  coprime  transfer  matrices  N,M  such  that  D  =  NM~l  (implicitly  we  have  that  M  is 
square  and  invertible).  The  following  theorem  addresses  the  existence  of  such  a  factorization. 

Theorem  1  Suppose  we  have  a  transfer  matrix  D  £  TZM,  and  (without  loss  of  generality)  that  a 
minimal  realization  for  D  is  given  by 


D  = 


A 

B' 

C 

D 

Ai 

0 

B\  ' 

0 

A-2 

b2 

,C\ 

c2 

D 

where  Ax  is  marginally  stable  and  A2  is  antistable.  Now  defining 


we  have  that  H  £  dom(Ric),  and  so  we  can  define 


X 

=  Ric(H)  <  0 

A-i 

BxBt2X 

Bi  1 

N  = 

0 

A2  +  B2BfX 

lb 

L  c'1 

C2  +  DB(X 

D  J 

(8) 


(9) 


(10) 

(11) 


5 


\A2  +  B2BfX 

b2] 

BjX 

I  _ 

(12) 


Then  we  have  that  D  =  NM  1 ,  with  N  marginally  stable,  M  inner,  and  (12)  a  minimal  realization 
of  M.  Furthermore  i/De  IZCoo  then  N  £  'R'H(y0,  and  N,M  are  right  coprime  over  IZHoo- 


Proof:  Note  first  that  H  clearly  has  no  imaginary  axis  eigenvalues  since  A2  is  antistable.  Now 
our  assumption  of  a  minimal  realization  for  D  implies  that  (A2,  B->)  is  controllable,  and  hence 
(A2,  B2Bf)  is  controllable.  Since  B2B2  is  positive  semidefinite,  these  facts  imply  that  H  £ 
dom(Ric),  so  X  =  Ric(H)  is  well  defined,  and  is  a  real  symmetric  matrix  satisfying  the  Riccati 
equation 

A^X  +  XA2  +  XB2B%X  =  0  (13) 

with  A2  +  B2B2  X  stable  (see  section  7.2  in  [Fra87]).  Rewriting  this  we  obtain 


(A2  +  B2Bt2  XfX  +  X(A2  +  B2Bt2  X)  -  XB2BT2X  =  0 


This  is  a  Lyapunov  equation,  and  so  the  fact  that  A2+B2Bf  X  is  stable  and  we  have  —XB2BfX  <  0 
implies  that  X  <  0.  By  construction  we  have  that  N  is  marginally  stable,  and  M  £  IZTtoo  is  square 
and  invertible.  Now  we  have 


\AX 

0 

B\BfX 

B1- 

0 

a2 

B2BfX 

b2 

0 

0 

a2  +  B2Bt2X 

B 2 

LCx 

c2 

DBlpT 

D 

Applying  the  state  transformation 


we  obtain 


[Ax 

BiBfX 

B\BfX 

Bi 

0 

a2  +  b2bix 

B2Bt2X 

b2 

0 

0 

a2 

0 

Lex 

C2  +  DBfX 

DBfX 

D 

and  hence  D  =  AM-1.  Furthermore  note  that 


[■Ax 

b,bt2x 

Bil 

0 

a2  +  b2b?x 

b2 

[Cl 

C2  +  DB'i'X 

D 

N 


M~M  = 

r  -A*  -  XB2BT2 

0 

-XB2BfX 
A2  +  B2B2  X 

-xb2- 

b2 

L 

BIX  | 

I 

_ 

'  -A\  -  XBoBf 

0 

at2x  +  xa2  +  XB2Bt2X 
A2  a  B2B?X 

0 

b2 

L 

0 

I 

' -A\  -  XB2Bl 

0 

0  ' 

0 

A2  +  B2Bf  X 

B2 

L 

0 

I J 

(where  we  made  use  of  the  state  transformation  T  =  y  ^  j  J  and  the  Riccati  equation  (13))  so 
that  M  is  inner.  Now  we  know  that  ( A2,B2 )  is  controllable,  A2  is  antistable,  and  A2  +  B2Bf  X 
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is  stable.  Applying  the  PBH  tests  to  (12),  and  using  these  three  facts,  it  is  easy  to  show  that 
( A2  +  B2B2  X,  B2)  is  controllable  and  ( B if X,  A2  +  B2B2  X)  is  observable,  so  that  (12)  is  a  minimal 
realization  of  M.  Finally  we  note  that  if  D  G  TZCoo  then  Ax  is  stable  and  hence  N  G  TZH^.  The 
fact  that  N,  M  are  then  right  coprime  over  TZTioo  may  be  verified  by  choosing  L  such  that  A  +  LC 
is  stable  (which  we  can  do  since  (C,  A)  is  observable)  and  checking  that 


P  = 


A  +  LC 

-B-LD  L 

0  B£X 

I  0  . 

G  'JZ‘H00 


(14) 


is  a  left  inverse  of 


□ 


This  theorem  provides  us  with  State  Space  formulae  to  construct  a  particular  factorization,  namely 
one  where  the  denominator  matrix  M  is  inner.  When  D  G  7 this  is  a  right  coprime  factorization, 
and  in  fact  it  is  easy  to  show  that  such  a  factorization  exists  if  and  only  if  D  G  TZCoo.  Note  from 
the  proof  that  for  D  £  TZCoo  Ax  has  imaginary  axis  eigenvalues  and  so  N  £  TZTToo  •  However  we  still 
have  an  inner  M,  and  a  marginally  stable  N,  such  that  D  =  NM~l,  with  P  G  TZT (see  (14))  a 

left  inverse  of 

The  existence  of  a  right  coprime  factorization  for  D  G  TZCoo  is  well  known.  What  is  new 
here  is  that  the  factorization  in  theorem  1  may  be  carried  out  for  arbitrary  D  G  TZM  and,  more 
importantly,  that  we  only  solve  a  Riccati  equation  of  the  same  size  as  A2,  i.e. ,  the  number  of 
antistable  poles  in  D.  Thus  we  only  invert  out  the  antistable  part  of  D  and  hence  obtain  a  minimal 
realization  for  M.  This  offers  clear  numerical  advantages  over  implementing  the  standard  formulae 
for  coprime  factorization  (see  [ZDG]  for  example),  which  would  require  solving  a  Riccati  equation 
of  the  same  dimension  as  A,  and  would  result  in  unobservable  modes  in  the  realization  for  M.  The 
numerical  advantages  of  the  formulation  presented  here  are  particularly  important  when  there  is 
pole  symmetry  about  the  imaginary  axis,  which  occurs,  for  example,  when  D  is  self-adjoint.  In 
that  case  the  Hamiltonian  matrix  from  the  standard  formulae  (see  [ZDG,  Fra87]  for  example) 


has  every  eigenvalue  repeated.  A  standard  approach  to  solving  the  Riccati  equation  relies  on 
computing  the  eigenspaces  of  the  Hamiltonian  matrix  H.  It  is  well  known  that  a  non-trivial  Jordan 
form  can  present  severe  problems  for  numerical  computation  of  the  eigenvalue  decomposition,  so 
that  a  straightforward  implementation  of  the  standard  formulae  for  coprime  factorization  could  run 
into  numerical  difficulties  on  this  type  of  problem. 

The  analogous  definitions  and  results  for  left  coprime  factorizations  follow  by  duality. 

Definition  5  Suppose  we  have  two  transfer  matrices  N,M  G  TZM..  Then  N,M  are  said  to  be  left 
coprime  (over  TZTtoa)  if  N,M  G  TZTTVJ,  they  have  the  same  number  of  rows,  and  ^  M  N  ^  is  right 
invertible  in  TZTioo- 

Given  a  transfer  matrix  D  G  TZM,  then  a  left  coprime  factorization  of  D  (if  it  exists)  consists  of 
two  left  coprime  transfer  matrices  N,M  such  that  D  = 
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Theorem  2  Suppose  we  have  a  transfer  matrix  D  G  1ZM,  and  (without  loss  of  generality)  that  a 
minimal  realization  for  D  is  given  by 


D  = 


'  A 

B  ' 

C 

D 

\Ai 

0 

Bi' 

0 

a2 

b2 

[c\ 

c2 

D 

where  Ai  is  marginally  stable  and  A2  is  antistable.  Now  defining 


(  Af  CJC2  \ 
V  o  -A2  ) 


(15) 


(16) 


we  have  that  H  G  dom(Ric),  and  so  we  can  define 

X  =  Ric(H)  <  0 


(17) 


Ai 

0 

1 

Bx 

xcfc. 

a2  +  XCJC2 

B2  +  XCfD 

Cx 

c2 

1 

D 

M  = 

r  a2 + xcjc2 

XCf] 

L 

I 

(18) 

(19) 


Then  we  have  that  D  =  M  1N ,  with  N  marginally  stable,  M  inner,  and  (19)  a  minimal  realization 
of  M.  Furthermore  if  D  G  IZCoo  then  N  G  IZHoo,  and  N,M  are  left  coprime  over  IZHoo ■ 


Proof:  Apply  theorem  1  to  DT  and  the  results  follow.  □ 


3.3  Inner-Outer  Factorization 

In  this  section  we  will  make  use  of  the  following  notions.  The  term  minimum  phase  refers  to  a 
system  (matrix)  with  all  its  zeros  (eigenvalues)  in  the  open  left  half-plane,  and  marginally  minimum 
phase  refers  to  a  system  (matrix)  with  all  its  zeros  (eigenvalues)  in  the  closed  left  half- plane.  The 
term  anti  minimum  phase  refers  to  a  system  (matrix)  with  all  its  zeros  (eigenvalues)  in  the  open 
right  half-plane. 

We  now  wish  to  develop  a  factorization  which  will  factor  out  the  anti  minimum  phase  part  of 
a  transfer  matrix.  It  is  useful  to  consider  another  special  type  of  transfer  matrix. 

Definition  6  Suppose  we  have  a  transfer  matrix  D  G  1ZM.  Then  D  is  outer  if  D  G  IZHoo  and  D 
is  right  invertible  in  IZHoo  ■ 

It  is  easy  to  show  that  a  square  transfer  matrix  D  G  1ZM  is  outer  if  and  only  if  D,D_1  G  IZHoo. 
Thus  an  outer  matrix  is  stable  minimum  phase,  i.e. ,  it  has  all  its  poles  and  zeros  in  the  open 
left  half-plane.  Note  that  in  an  analogous  fashion  we  can  define  a  matrix  D  G  1ZM  as  co-outer  if 
D  G  IZHoo  and  D  is  left  invertible  in  IZHoo-  It  is  clear  that  for  square  transfer  matrices  D  is  outer 
if  and  only  if  D  is  co-outer. 

If  we  can  factor  a  matrix  D  G  TZM.  as  D  =  DjD0,  with  D <  inner,  and  D<,  outer,  then  this 
is  referred  to  as  an  inner-outer  factorization  of  D.  Clearly  we  must  have  D  G  IZHoo  for  such  a 
factorization  to  exist,  and  in  fact  it  can  be  shown  that  such  a  factorization  exists  for  all  D  G  IZHoo 
(see  [ZDG,  Fra87]).  Here  we  need  only  concern  ourselves  with  a  special  case  of  this  factorization, 
namely  when  D  is  square  and  invertible. 


Theorem  3  Suppose  D  £  TZM  is  square,  invertible,  and  marginally  stable.  Further  suppose  (with¬ 
out  loss  of  generality)  that  we  have  a  minimal  realization  for  D-1 


D 


1  _ 

A 

B' 

_C 

D 

A\ 

0 

Bb 

0 

A2 

b2 

.Cl 

C2 

D 

with  A\  marginally  stable,  and  A2  anti-stable,  so  that  defining 

A  =  A  —  BD~lC 
B  =  BD~l 
C  =  ~D-lC 

D  =  d-1 


(20) 


(21) 


we  have  that 


A 

B 

A 

B 

c 

D 

4  c2 

D 

is  a  minimal  realization  for  D.  Now  defining 


we  have  that  H  £  dom(Ric),  and  so  we  can  define 


(22) 


(23) 


X  =  Ric(H)  <  0 


(24) 


r  a2  +  B2if(x 

B2 

b2  X 

I  _ 

A 

b' 

Cl  C2  —  BTf  X 

D 

(25) 

(26) 


Then  we  have  that  D  =  DjD0  with  D*  inner,  and  (25)  a  minimal  realization  for  D;,  and  D0 
marginally  stable,  marginally  minimum  phase.  Furthermore  if  D  £  RTtoo  then  D0  £  TTHoo  and  if 
D”1  £  RCoo  then  De  is  (right)  invertible  in  7 ZH^. 


Proof:  Apply  theorem  1  to  factor  D  1  as  D  1  =  NM  1  with 


rAi 

BxB\X 

Bx  1 

A  = 

0 

A2  +  B2BJ  X 

b2 

[Cl 

C\  +  DBfX 

D  J 

(27) 


r  A2  +  B2BfX 

b2] 

[  BfX 

I  _ 

(28) 


where  X  =  Ric(H)  <  0.  Note  that  N  is  marginally  stable,  M  is  inner,  and  (28)  is  a  minimal 
realization  of  M.  Note  also  that  N  is  square  and  invertible  (as  is  M).  Now  define  D,-  =  M  and 
D0  =  A-1,  and  we  have 


D  =  (D-1)-1  =  (AM-1)-1  =  MN~l  =  D,D0 
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with  D,  inner,  (25)  a  minimal  realization  of  D,-,  and  D0  marginally  minimum  phase  (since  D0  1  is 
marginally  stable).  Now  by  construction  we  have 


ByB'fX 

B il 

Dc  = 

0 

A2  +  B2BlX 

b2 

Lei 

C2  +  DBfX 

D  J 

A 

B  ' 

.  C\  C2-  BIX 

D 

(29) 


so  that  D0  is  marginally  stable.  Finally  note  that  if  D  G  TZHoo  then  A  is  stable  so  D0  €  TZTioo,  and 
if  D_1  G  TIjCoo  then  N  G  TZTioo,  so  D0  is  (right)  invertible  in  TZTioo  (with  D"1  =  N).  □ 


In  the  case  that  D  G  TZTioo,®*1  G  TZCoo  we  have  that  D0  is  outer  and  hence  this  is  exactly 
an  inner-outer  factorization.  Note  once  again  that  although  the  inner-outer  factorization  is  well 
known,  we  are  here  only  solving  a  Riccati  equation  of  the  same  dimension  as  Ay ,  i.e. ,  the  anti 
minimum  phase  zeros  of  D.  Thus  by  only  factoring  out  the  anti  minimum  phase  part  of  D  we 
obtain  a  minimal  realization  for  D*.  This  offers  numerical  advantages  over  the  standard  approach, 
and  similar  comments  to  those  made  in  section  3.2  apply. 

Once  again  we  may  define  a  dual  version  of  this  factorization,  namely  if  we  can  factor  a  matrix 
D  G  TZM  as  D  =  DcoDci,  with  Dci  co-inner,  and  Dco  co-outer,  then  this  is  referred  to  as  a  co- 
outer-co-inner  factorization  of  D.  Again  the  results  follow  easily  by  duality,  and  we  do  not  include 
them  here 


3.4  Combined  Factorizations 

The  machinery  presented  in  the  preceding  subsections  enables  us  to  carry  out  a  combined  factor¬ 
ization,  which  splits  off  the  anti-stable  poles  and  anti  minimum  phase  zeros  of  D  using  all  pass 
functions. 


Theorem  4  Suppose  we  have  a  square  invertible  transfer  matrix  D  G  TZM. 
factored  as 


D  —  DapaDapi  ®  smp 


Then  D  may  be 


(30) 


where  Dapa  is  anti  stable  all  pass,  Daps  is  stable  all  pass ,  and  Dsmp  is  marginally  stable,  marginally 
minimum  phase.  Furthermore  if  D  G  TZC^  then  Dsmp  is  stable,  and  if  D-1  G  TZCoo  then  Dsmp  is 
minimum  phase. 


Proof:  Apply  theorem  2  to  obtain 

D  =  M~1N 

with  M  inner,  and  N  marginally  stable.  Define  Dapa  =  M~l  and  it  is  easy  to  check  from  (19)  that 
M~l  is  anti-stable.  Now  apply  theorem  3  to  obtain 


N  —  n  D 

ap$  smp 


with  Daps  inner  and  Dsmp  marginally  stable,  marginally  minimum  phase.  The  results  now  follow. 

□ 


Note  that  for  D,D  1  G  TZCoo  (which  will  be  the  case  when  this  result  is  employed  in  the  mixed  p 
synthesis  problem)  this  allows  us  to  factor  D  such  that 

D  =  DapDsmp  (31) 
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where  Dap  =  DapaDapil  is  an  all  pass  function,  and  Dsmp  is  stable  minimum  phase.  The  State 
Space  formulae  to  perform  this  factorization  are  given  in  theorems  2  and  3.  In  order  to  compute 
this  factorization  we  need  to  solve  two  Riccati  equations.  Note  that  these  two  Riccati  equations 
have  dimension  equal  to  the  number  of  anti-stable  poles  of  D,  and  the  number  of  anti  minimum 
phase  zeros  of  D  respectively. 


3.5  Spectral  Factorization 

Finally  we  consider  spectral  factorization,  which  is  essentially  a  transfer  matrix  equivalent  of  the 
square  root.  We  are  concerned  with  a  particular  form  of  spectral  factorization,  that  arises  from  the 
mixed  p,  upper  bound. 


Theorem  5  Suppose  we  have  a  square  transfer  matrix  G  G  TIM..  Further  suppose  (without  loss  of 
generality)  that  a  minimal  realization  for  G  is  given  by 


A 

B' 

C 

D 

Now  define  the  symmetric  positive  definite  matrices  R  and  Q  as 


R  =  I+DTD 
Q  =  I  +  DDT 


and  define  the  Hamiltonian  matrix 

it  _  (  A  -  BR'1DtC  -BR~lBT 
V  -CtQ-xC  -{A- BR~1DtC)t 

Then  H  G  dom(Ric)  so  we  may  define 


X  =  Ric(H)  >  0 


'  A  -  BR~1DTC  -  BR~1BTX 

BR~i  ' 

^  h  — 

-R~\DTC  +  BTX) 

R-> 

Then  we  have  that  Gh  €  TTH^  and 

(I+G~G)~l  =GhG~ 

with  GGh  G  IZHoo  given  by 

II 

'  A  -  BR~1DTC  -  BR~lBTX 

BR~ i 

Q~1(C  -  DBT X) 

DR-i 

Proof:  First  we  make  the  definitions 

A  =  a-br~1dtc 

B  =  BR~ i 

C  =  Q~lC 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 
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It  is  easy  to  show  that  {A,  B)  is  controllable  and  (C,  A)  is  observable.  But  now  we  note  that  II 
has  the  form 

H  = 

and  hence  we  have  that  H  G  dom(Ric),  so  that  X  =  Ric(H)  is  well  defined  and  X  >  0  (see  lemma 
3  in  [DGKF89])  satisfies  the  Riccati  equation 

(A  -  BR~1DtC)tX  +  X(A  -  BR~1DtC )  -  XBR~1BTX  +  CTQ~1C  =  0  (39) 


with  A- BR~lDTC  - 
invertible.  Now  define 


BR  1BT X  stable.  Thus  by  construction  Gh  G  RHoo, 


A 

B 

[r-Hdtc  +  BTX) 

Rl*  . 

and  Gh  is  square  and 


(40) 


and  we  have 


N~N  = 

A 

—(CTD  +  XB)R~1(DTC  +  BTX) 

o  ^ 

1 

B 

-CTD  -  XB 

[  DTC  +  IVX 

B1 

R  \ 

A 

0 

B 

-( ATX  +  XA  +  CTC ) 

-AT 

—CTD  —  XB 

[  DTC  +  B‘X 

B'r 

R  \ 

'  A 

0 

B  ‘ 

-CTC 

-AT 

-CTD 

[  D'C 

Br 

R  J 

I  +  G~G 


Note  we  have  made  use  of  the  Riccati  equation  (39)  and  the  similarity  transformation  T  = 


/ 

-X 


Thus  we  have 


(I  +  G~G)~X  =  ( N~N)~l  =  N-'N-1-  =  GhGh 


Finally  we  have  that 


'  A  -  BR~1DTC  -  BR~1BTX 

0 

BR-i  ' 

GGh  — 

-BR~1DtC  -  BR~1BTX 

A 

BR~i 

-DR~l(DTC  +  BTX) 

C 

DR ~i 

A  0 

0  ' 

—BR~lDTC  -  BR~1BTX  A  -  BR~lDTC  -  BRrlBTX 

BR~i 

L  -DR-1(DTC  +  BTX)  C  —  DR~l{DT  C  +  BTX) 

DR- it  J 

'  A  -  BR~1DTC  -  BR~lBTX 

BR~i  ' 

Q~l(C  —  DBT  X) 

DR~i 

(where  we  have  used  the  similarity  transformation  T 


j  ^ )  and  hence  GGh  G  RHoo  ■ 


□ 
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Essentially  Gh  is  a  stable  realization  of  (I  +  G~G)~%.  Note  that  we  do  not  make  any  restrictions 
on  the  pole/zero  locations  of  G  for  the  existence  of  this  factorization.  In  particular  G  may  have 
poles  and/or  zeros  on  the  imaginary  axis,  and  we  still  obtain  stable  realizations  for  Gh  and  G G h • 
It  will  be  seen  later  that  this  degree  of  generality  is  required  for  the  mixed  fi  synthesis  problem. 

4  Controller  Synthesis 

We  now  wish  to  consider  the  problem  of  synthesising  a  controller  to  meet  a  certain  robust  perfor¬ 
mance  requirement.  First  we  place  the  problem  in  a  standard  framework  as  shown  in  figure  1. 


Figure  1:  Feedback  interconnection  for  /i  synthesis 

The  transfer  matrix  P  is  assumed  here  to  contain  all  the  appropriate  weights  for  the  problem,  as 
well  as  the  interconnection  structure  pertaining  to  the  model  of  the  physical  plant.  Note  that  we 
can  define  the  transfer  matrix  M(P,  K),  by  absorbing  the  controller  K  into  P  as  shown  in  figure  1, 
and  so  place  the  problem  in  the  standard  robust  performance  //,  analysis  framework  (see  [PD93]). 

We  are  interested  in  developing  methods  to  tackle  the  general  /i  synthesis  problem,  namely  that 
of  finding  a  controller  K  achieving 

inf  sup  /i*:(M(P ,  K)(jw))  (41) 

Wfe/t-s  uiSIR 

where  K,s  denotes  all  real  rational  proper  controllers  that  nominally  stabilize  P  (i.  e. ,  render 
M(P,K)  internally  stable).  The  control  engineering  motivation  for  this  problem  is  extensive,  and 
may  be  found  elsewhere  (see  [PD93,  YND91,  You93]  and  the  references  therein),  together  with 
the  precise  mathematical  robustness  interpretation  of  (41).  Here  we  merely  note  that  a  solution 
to  (41)  gives  us  a  controller  achieving  the  best  possible  robust  performance  (in  an  7i.rX}  norm 
sense)  with  respect  to  a  certain  perturbation  class.  Note  also  that  here  we  are  explicitly  including 
real  uncertainties  in  the  block  structure,  so  that  this  synthesis  procedure  considers  parametric 
uncertainty  directly,  rather  than  simply  approximating  the  real  uncertainties  with  complex  ones, 
and  applying  standard  complex  fi  synthesis  techniques. 
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4.1  Complex  p  Synthesis  and  “ D-K  Iteration” 

Prior  to  tackling  the  mixed  p  synthesis  problem,  we  first  include  a  very  brief  review  of  complex  p 
synthesis  (where  mr  —  0).  For  a  more  detailed  exposition  see  [Doy85,  SD91].  First  note  that  the 
synthesis  problem  (41)  is  not  tractable,  and  so  we  consider  the  problem  given  by  replacing  complex 
p  by  its  upper  bound  in  (41),  namely  (see  [Doy85]) 

inf  sip  inf  <r(B(u>)M(P,K)(j. o)D~\u))  (42) 

KeXs  a/g®  d(w)£T>k. 

Note  that  if  we  fix  K  then  the  problem  of  finding  D(u )  is  just  the  standard  complex  p  upper  bound 
problem  (across  frequency)  which  is  a  convex  problem  and  can  be  efficiently  solved.  If  we  choose 
D(u> )  matrices  at  a  set  of  frequency  points  (from  the  p  upper  bound)  we  can  fit  a  real  rational, 
stable,  minimum-phase  transfer  matrix  to  them.  If  we  fix  this  D  transfer  matrix  then  the  problem 
of  finding  K  reduces  to  a  standard  'H<x>  problem,  as  follows.  Consider  the  feedback  interconnection 


Figure  2:  Standard  framework  for  Hoo  optimal  control 

in  figure  2.  Assume  once  again  that  we  have  collected  all  the  relevant  quantities  in  the  nominal 
generalized  plant  P,  and  denote  by  Tzw  the  closed  loop  transfer  function  from  w  to  z.  Then  the 
problem  of  choosing  a  real  rational  proper  controller  K,  so  as  to  render  this  feedback  interconnection 
internally  stable,  and  minimize  UT^Uoo,  is  the  standard  Hoo  optimal  control  problem.  This  problem 
is  convex  and  a  solution  was  obtained  in  [DGKF89],  for  which  commercial  software  is  now  available 
[BPD+91]. 

The  above  approach  leads  to  the  following  “ D-K  iteration”  scheme,  which  attempts  to  find  a 
complex  p  optimal  controller  (and  hence  the  best  achievable  robust  performance  for  problems  with 
structured  complex  uncertainty): 

Procedure  1  ( D-I(  Iteration) 

1.  Find  an  initial  estimate  of  the  scaling  matrices  D(uj)  pointwise  across  frequency.  One  possi¬ 
bility  is  to  use  the  identity  matrix  at  each  point. 

2.  Find  a  State  Space  realization,  D,  fitting  the  pointwise  scaling  matrices  D(u)  with  a  stable 
minimum  phase  system  (so  that  D  and  D-1  are  stable).  Augment  this  with  an  identity 
matrix  of  the  appropriate  size  so  that  D  is  compatible  with  P.  Construct  the  State  Space 
system  Pd  =  DPD-1. 
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3.  Find  the  optimal  controller  K  minimizing  ||M(Pd,  K)||oo  over  all  all  stabilizing,  proper, 
real  rational  controllers  K. 

4-  Find  D(u)  solving  the  minimization  problem 

inf  o(D{u)M(Y>,K)(iu)D-\u)) 

D(w)et>ic 

pointwise  across  frequency. 

5.  Compare  D{u)  with  the  previous  estimate  D(u>).  Stop  if  they  are  close,  else  replace  D(u) 
with  D(u>)  and  return  to  step  2. 

This  iteration  (assuming  perfect  State  Space  realizations  of  D{ u>))  is  monotonically  nonincreasing, 
so  that  we  are  guaranteed  convergence  to  a  local,  minimum  of  the  problem.  Having  converged  to 
such  a  point  the  controller  K  from  step  3  is  the  resulting  complex  p  synthesis  controller.  Note  that 
although  the  individual  problems  (the  //,  upper  bound  and  7icx>  optimal  control)  are  convex,  the 
joint  problem  (/1  optimal  control)  is  not  convex  (see  [Doy85]  for  a  counterexample  to  convexity). 
Thus  the  UD-K  iteration”  described  above  is  not  guaranteed  to  converge  to  the  global  optimum  of 
the  problem.  However  many  designs  have  been  performed  using  this  technique  in  recent  years  (see 
[BCD89,  BD90]  for  example),  and  it  has  usually  been  found  to  work  well  in  practice.  There  are 
further  subtleties  to  the  above  procedure,  which  we  will  not  go  into  here  (see  [SD91,  BDG+91]  for 
more  details). 

This  “ D-K  iteration”  scheme  is  the  basis  of  complex  p  synthesis  controller  design.  The  design 
process  consists  essentially  of  the  following  stages. 

Procedure  2  (Complex  p  Synthesis) 

1.  Decide  on  an  appropriate  interconnection  structure  to  model  the  system,  including  the  uncer¬ 
tainty  structure  against  which  robustness  is  desired. 

2.  Choose  appropriate  weights  to  reflect  the  desired  performance  specifications,  and  any  infor¬ 
mation  known  about  the  uncertainties. 

3.  Implement  the  above  “D-K  iteration”.  Note  that  this  involves  deciding  on  an  appropriate 
frequency  range  of  interest,  and  selecting  the  order  of  the  State  Space  fit  in  step  2. 

This  procedure  is  by  no  means  mechanical,  and  a  great  deal  of  engineering  judgement  is  still 
required,  particularly  in  steps  1  and  2.  Also,  as  with  any  control  design  technique,  the  full  design 
process  may  well  involve  an  iterative  application  of  the  above  process,  together  with  an  evaluation 
of  each  resulting  controller. 


4.2  Mixed  p  Synthesis 

As  in  the  complex  case  we  note  that  (41)  is  not  tractable,  and  so  we  consider  the  problem  given 
by  replacing  p  by  its  upper  bound  in  (41).  Now  for  any  matrix  M  G  C”xn  We  have  that  (see 
[FTD91,  You93]) 


pAM)  <  inf  inf  <  8  :  0 

V  '  —  DeT>K  ,G£Qk.  0€K,/3>O  I 


/  DMD~l 

V  p 


-  jG  )  (I  +  G2)A  \  <  1 


(43) 
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Note  that  if  we  enforce  the  choice  G  =  0„  then  we  recover  the  standard  complex  p  upper  bound. 
The  additional  degrees  of  freedom  we  have  in  choosing  the  G  matrix  enable  us  to  exploit  the  phase 
information  we  have  about  the  real  uncertainties  to  get  a  better  bound.  Note  also  that  this  problem 
may  be  rearranged  into  a  quasi-convex  optimization  problem,  for  which  efficient  computational 
software  exists  (we  refer  the  interested  reader  to  [YND92]  for  details  of  the  practical  computation 
of  the  mixed  p  upper  bound  and  associated  scaling  matrices).  This  leads  us  to  consider  the  following 
synthesis  problem. 


inf 
k  exs 


sup  inf  inf  {/3(«)  :  T  <  1} 

we®  %)6®ic,<s(u)eeic  ^(w)e®,/?(w)> o^v 


(44) 


where 


VV  /?(«) 


-jG(w))  (J  +  Ga(w))~4) 


(45) 


The  basic  idea  behind  our  approach  to  mixed  p  synthesis  is  evident  from  these  equations  (44,45). 
For  fixed  K  then  the  problem  of  finding  D(u>),  G(u>),(3(u)  is  just  the  mixed  p  upper  bound  prob¬ 
lem  (across  frequency),  which  is  a  convex  problem  and  can  be  efficiently  solved.  Having  chosen 
D(u),G(tjj),  f3{u>)  at  a  set  of  frequency  points  (from  the  p  upper  bound),  we  will  fit  real  rational 
transfer  matrices  to  them,  in  such  a  way  that  the  interconnection  is  stable.  For  fixed  D,G,/? 
transfer  matrices  then  the  problem  of  finding  K  will  be  reduced  to  a  standard  Moo  problem  which 
again  is  convex  and  can  be  efficiently  solved.  This  will  lead  to  a  “ D,G-K  iteration”  scheme,  for  the 
mixed  p  synthesis  problem. 

Of  course  we  will  suffer  from  the  same  limitations  as  the  “ D-K  iteration”  procedure  for  complex 
p  synthesis,  namely  that  the  joint  problem  of  optimizing  D,G  and  K  is  not  convex.  Thus  the 
scheme  will  not  guarantee  to  find  the  global  optimum  of  the  problem,  but  rather  only  a  local 
minimum.  Nevertheless  we  will  obtain  a  procedure  which  exploits  the  phase  information  in  the 
real  uncertainties,  via  the  G  scaling  matrix,  and  so  has  the  potential  of  greatly  improving  on  the 
present  approach,  where  one  simply  covers  the  real  uncertainties  with  complex  ones,  and  applies 
complex  p  synthesis. 


4.2.1  UD,G-K  Iteration” 

Note  that  the  problem  statement  in  equations  (44,45)  is  not  quite  in  the  form  of  a  (scaled)  singular 
value  minimization,  but  rather  a  minimization  subject  to  a  (scaled)  singular  value  constraint.  In 
order  to  proceed  we  need  to  exploit  some  properties  of  the  mixed  /. i  upper  bound. 

Theorem  6  Suppose  we  have  matrices  M  G  Cixn  ,D  G  VK,G  G  Gk.  and  a  real  scalar  f3  >  0  such 
that 

a)  V +&)-*)<! 

Then  for  any  real  (3  >  f3  there  exist  matrices  D  G  T>k  and  G  G  Gk  smc/j  that 


a 


^ DMD -1 


■j  G  (i  +  G2y 


<  1 
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Proof:  For  convenience  define  M  =  DMD  *,  and  we  have 

a((f -jG](/  +  G2)^)  <  1 

—  Of+JS)  (f  -&)  <  I  +  G2 

— >  M*M  +j(3(GM  -  M*G)  <  (32I  <  /32/ 

—  ^+j  $(GM-M*G)  <  I 

Thus  defining  D  ±  D,  G  =  we  have 

M  +  X(GM  -  M*G)  <  I 

—  (x+j^)  (f -j^)  <  I  +  d2 

— ►  o[{^fX-iG){I  +  G2Y^)  <  1  □ 

Theorem  7  Suppose  we  have  matrices  M  G  C 7lxn,D  G  Pjt,G  G  nnd  reed  scalars  f3  >  0  and 

0  <  r  <  1  such  that 

KFF-jG)(/+Gn>sr 

Then  there  exists  D  G  V>c,G  G  (7/e  suc/i  that 


Proof:  Now  define  M  =  DM^  and  we  have 

a((M-jG')(/  +  G2)-i)  <  r 

— ►  (M*+jG)(M-jG)  <  r2(/  +  G2) 

— M*M  +i(GM  -  M*G)  <  r2I+(r2-l)G2  <  r2I 
— »■  ^  +  <  / 

so  that  defining  D  =  D,  G  =  ®  we  have 

^  +  X(GM  -  M*G)  <  I 

—  fe+jG)  (f  -jG)  <  ^  +  G2 

— >  a((^f^-jG')(/  +  G2)-^)  <  1  □ 

These  two  theorems  tell  us  something  about  how  the  mixed  p  upper  bound  scales.  Note  that 
both  of  these  results  are  trivial  in  the  complex  case,  where  G  =  0„.  For  the  mixed  case  however 
the  results  are  not  so  obvious,  and  in  particular  it  is  important  to  note  that  theorem  7  applies  only 
if  r  <  1.  From  this  theorem  we  see  that  whenever  we  have  the  appropriate  scaling  matrices  and 
r  <  1,  then  r/3  is  an  upper  bound  for  pK(M).  However  when  r  >  1  we  cannot  conclude  that  r/3  is  an 
upper  bound  in  the  mixed  case,  even  though  this  conclusion  is  obviously  true  in  the  complex  case. 
In  fact  it  is  easy  to  construct  mixed  examples  with  r  >  1  where  r/3  is  strictly  less  than  pYM)- 
We  now  propose  the  following  “ D,G-K  iteration”  for  the  mixed  p  synthesis  problem: 
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Procedure  3  ( D,G-K  Iteration) 

1.  Find  initial  estimates  of  the  scaling  matrices  D{u),  G(oj),  and  the  real  positive  scalar  (3*.  One 
possibility  for  the  scalings  D(us)  is  to  choose  them  as  the  identity  matrix  at  each  frequency 
point.  If  G(u)  is  chosen  to  be  the  zero  matrix  at  each  frequency  point,  then  /3*  is  arbitrary 
(so  choose  say  /?,  =  1),  otherwise  it  must  satisfy 

a  ( _  jG(o>))  ( I  +  G(o;)2)-^  <  1 
for  all  oj,  for  some  stabilizing  controller  K. 

2.  Fit  State  Space  realizations  D  and  G  to  the  pointwise  scaling  matrices  D(u)  and  }G(u>),  so 
that  D(jw)  approximates  D{lo),  and  G(jw)  approximates  ]G{u>).  Now  replace  D  and  G  with 
appropriate  factors  so  that  D,  D-1,  G&,  and  GGh  are  all  stable,  where  Gh  is  a  spectral  factor 
satisfying  (/ +  G~G)“1  =  GaGJ“.  Augment  D  and  Gh  with  identity  matrices,  and  G  with 
a  zero  matrix,  of  appropriate  dimensions  so  that  D,  G,  Gh  are  all  compatible  with  P.  Form 
the  State  Space  system  PDG  =  (DPD-1  —  /3*G)Gft. 

3.  Find  the  optimal  controller  K  minimizing  ||M(PDG,  K)||oo  over  all  all  stabilizing,  proper, 
real  rational  controllers  K. 


4-  Compute  /?*  as 


/?*  =  sup  inf  inf  :  T  <  1} 

wgK  D(u)eT>ic,GOeeic  P(u)eB.,f}(w)>0 


where 


T  =  a 


^  D(u>)M(P,K)Qfj)D-X<j) 


P(u>) 


iG{u)\(I+G\u))~ 


5.  Find  D(u),G(lv)  solving  the  minimization  problem 


inf  a 
D(w)eVK,G(u)eGic 


,//^(«)M(P,K)au;)i)-V) 


/?* 


-jG(u)  (7  +  G2H)-i 


pointwise  across  frequency. 

6.  Compare  the  new  scaling  matrices  D{u),G(ui)  with  the  previous  estimates  D(u>),  G(lo).  Stop 
if  they  are  close,  else  replace  D(uj),  G(co)  with  D(l>),G(u)  respectively  and  return  to  step  2. 

We  will  see  that  this  iteration  (assuming  perfect  State  Space  realizations  of  D(u)  and  (^(w))  is 
monotonically  nonincreasing,  so  that  we  are  guaranteed  convergence  to  a  local  minimum  of  the 
problem.  Having  converged  to  such  a  point  then  the  controller  K  from  step  3  is  the  resulting  mixed 
ji  synthesis  controller,  and  it  satisfies 

sup  /r/c(M(P ,  K)(jw))  <  /?* 

Note  that,  as  in  the  complex  p  case,  we  use  the  Hoo  optimal  control  solution  to  synthesize  the 
controller.  By  construction  the  scaled  system  PDG  fits  into  the  standard  H optimal  control 
framework,  as  shown  in  figure  3. 
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Figure  3:  The  UK”  iteration  of  the  “ D,G-K  iteration”  procedure 


The  fact  that  this  iteration  is  monotonically  nonincreasing  follows  from  theorems  6  and  7  in 
the  following  way.  Suppose  we  have  run  through  several  cycles  of  the  above  procedure,  and  have 
arrived  at  step  3,  where  we  wish  to  design  a  new  controller.  Since  /?*  was  defined  in  (the  previous) 
step  4  as  the  maximum  across  frequency  of  the  upper  bound  we  have  that  /3,  >  /3(u)  for  all  u, 
with  =  f3(oj)  only  for  those  frequency  points  whose  upper  bound  achieves  the  maximum  across 
frequency.  Then  theorem  6  guarantees  us  that  in  (the  previous)  step  5  we  could  achieve 

v  (/  +  G’M)-*)  <  1 

for  all  frequency  points  with  /?»  >  /?(w)  and 

U  P* 

for  points  with  /?*  =  /?(o>).  Now  our  previous  controller  from  step  3,  K,  stabilized  the  inter¬ 
connection  M(P,K),  and  hence,  by  our  choice  of  factorizations  for  D  and  Gh,  K  stabilizes 
M(Pdg,K)  as  well.  Note  that  the  above  formulae  imply  that  K  achieved  ||M(PDG,  K)^  =  /?„, 
so  that  we  are  guaranteed  the  existence  of  a  stabilizing  controller  achieving  at  least  this  level  for 
||M(Pdg,Km)||m.  Thus  our  new  7foo  optimal  controller  in  step  3  will  satisfy  ||M(PDG,  Kn™)^  = 
r/3*  with  r  <  1.  But  theorem  7  implies  that  rfi.,  is  now  an  upper  bound  for  //.  across  frequency,  and 
so  the  new  computation  of  (3*new  in  step  4  will  yield  P*new  <  r/3 *  <  /?».  Thus  the  iteration  is  mono¬ 
tonically  nonincreasing,  and  furthermore  if  the  new  W^-optimal  controller  from  step  3  achieves 
any  improvement  over  the  old  one  (i.e. ,  r  <  1)  then  the  peak  value  of  the  fi  upper  bound  across 
frequency  is  strictly  decreased. 

Note  also  that  the  Hoc  norm  of  the  interconnection  M(PDG,K)  (formed  by  closing  the  new 
interconnection  PDG  from  step  2  with  the  previous  controller  K)  is  determined  by  those  frequencies 


jG(W)  (/+(?»)- 
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which  have  the  maximum  value  for  the  mixed  n  upper  bound.  This  places  the  emphasis  for  the 
7foo  norm  minimization  at  the  frequencies  where  mixed  /i  is  largest  (which  is  just  what  we  want). 
Furthermore  the  way  we  have  constructed  the  interconnection  Pdg  guarantees  that  if  the  new 
controller  Kneu)  reduces  the  Moo  norm  of  M(PDG,K„ew),  then  the  peak  value  of  the  g,  upper 
bound  across  frequency  is  reduced. 

Note  that  in  step  4  we  compute  scalings  D(u),G(u>)  which  minimize  the  scalar  /3(lo)  at  each 
frequency  subject  to  a  (scaled)  singular  value  constraint,  and  this  is  exactly  the  mixed  fi  upper 
bound  computation  across  frequency.  A  closely  related  problem  is  solved  in  step  5,  where  we 
compute  scalings  D(u>),G(u )  to  minimize  a  (scaled)  singular  value  problem  with  /?*  fixed  across 
frequency.  If  we  consider  complex  problems  then  the  separate  computation  of  D(u),G(u)  in  step 
4,  and  then  D(u>),  G(u>)  in  step  5  is  unnecessary,  since  it  is  easy  to  see  that  the  two  sets  of  optimal 
scalings  may  be  taken  to  be  the  same.  For  mixed  problems  however  the  two  sets  of  scalings  in  steps 
4  and  5  can  be  quite  different,  and  it  is  important  to  use  the  scalings  D(w),G(u)  from  step  5  for 
fitting  (in  step  2). 

4.2.2  Forming  The  New  Interconnection 

In  order  to  implement  the  “ D,G-K  iteration”  described  in  the  preceding  subsection  we  need  to 
be  able  to  carry  out  certain  factorizations,  required  in  step  2.  These  factorizations  are  used  to 
ensure  that  if  K  stabilizes  M(P,K),  then  K  also  stabilizes  M((DPD_1  —  /?*G)Gft,K),  and  vica 
versa  (for  the  converse  direction  we  need  a  slight  technical  assumption  that  we  don’t  get  right 
half-plane  pole/zero  cancellations  when  forming  DPD_1G/,).  The  fact  that  K  stabilizes  both  these 
interconnections  was  used  when  showing  that  the  aD,G-K  iteration”  is  monotonically  nonincreasing 
(and  hence  converges). 

The  problem  of  fitting  pointwise  frequency  data  D(lo),G(u>)  with  State  Space  realizations  D,  G 
will  be  considered  in  the  next  subsection.  For  now  let  us  assume  that  we  can  obtain  such  fits.  We 
will  see  in  section  4.2.3  that  we  may  assume  that  D  is  square  and  invertible  with  D,D-1  E  TZCoo, 
but  the  only  assumption  we  will  make  about  G  is  that  it  is  square. 

Given  State  Space  realizations  for  D  and  G  we  want  to  factor  them  as  outlined  in  step  2  of 
the  “ D,G-K  iteration,”  described  in  procedure  3,  so  that  for  the  new  D,  G,  G^  the  interconnection 
M((DPD“1  —  /?* G)G/,,K)  is  stable,  and  we  would  like  to  do  this  without  affecting  the  value  of 

a(M((D(ju;)P(ju;)D-1(j^)  -  &G(ju,))Gft(ju,),  K(jW))) 

for  any  frequency.  This  can  be  accomplished  by  using  the  machinery  developed  in  section  3,  applied 
to  each  block  of  D  and  G. 

First  of  all  suppose  we  have  fitted  D(uj)  with  a  transfer  matrix  D,  but  have  not  yet  fitted  G{ui). 
Apply  theorem  4  to  each  block  of  D  to  factor  out  an  all  pass  function  so  that  the  resulting  block 
of  D  is  stable  minimum  phase  (hence  D  and  D ~ 1  are  stable).  Thus  we  have  a  new  transfer  matrix 
Dsmp  which  is  stable  minimum  phase,  and  at  each  frequency  point  we  have 

D(jw)  =  l7(w)D.m|,(jw)  (46) 

where  U(cj)  is  a  unitary  block  diagonal  matrix  (with  the  same  block  structure  as  D).  Then  denote 
D(jw),  U (w)  and  Dsmp(ju>)  by  Dw,  Uw  and  Dsw  respectively.  Finally  denote  the  pointwise  frequency 
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(47) 


data  G(w)  (to  be  fitted)  by  Gw,  and  we  have  the  following  equivalences: 

°  f^M(P,K)(jo;)i?-1  -  j/3tGw)  (/  +  G>)-i) 

=  aUuwDswM(P,K)(ju)D^U:  -  j/3*cO  (/  +  G2W )-*) 

=  a  f  f2?#tBM(P>K)(jw)JDrJ  -  j P*U*WGWUW)  U*(I  +  Gl)~^uJ) 

=  a  {{DswM(P,K)(jU)D^  -  j /?,<?«)  (I  +  <?>)-*) 

where  Gw  =  U*GWUW.  Note  that  Gw  £  Gk,  so  now  choose  the  transfer  matrix  G  so  that  G(jw) 
is  a  fit  to  j Gw  (rather  than  j Gw).  We  may  now  apply  theorem  5  to  each  non-zero  block  of  G  to 
construct  Gft,  with  Gfe  and  GGft  stable,  and  Gh  having  the  same  block  diagonal  structure  as  G. 
Furthermore  at  each  frequency  G /,  satisfies 

(7  +  G~(jn;)G(ju;))-1  =  {I  +  Gl)~l  =  Gft(jo;)G^(ja;)  (48) 

Now  for  convenience  denote  DSWM(P, K)(jw)7?jJ  by  MDw ,  and  Gft(ju>)  by  Ghw .  Then  note  that 
we  have 

<72  —  j/?*Gu>)  (. r  +  Gl)- 4) 

=  \((I+Gl)-i  (M*Dw-p*G*Qu))(MDw-/3*G(ju))(I  +  Gl)-i) 

=  J{(M*Dw  -I3>G*Qu))(MDw  -PM^))(I  +  Girl)  (49) 

=  \((M*Dw  -  /?*G*(ja>))  (MDw  -  /3*G(j u))Ghv>GZj 
=  X  (GZW  (M*Dw  -  P*G*Qu))  ( MDw  -  P.G(fr))  Ghw) 

—  o2  ((Mdw  —  P*G(juj))Ghw) 

so  that  this  factorization  does  just  what  we  require,  and  the  appropriate  State  Space  formulae  for 
D,  Gh  and  GGi  may  be  obtained  by  applying  the  results  from  section  3  to  each  block  of  D  and 

G. 

Note  that  the  final  step  of  this  process  is  to  augment  D  and  G/,  with  an  identity  matrix,  and  G 
with  a  zero  matrix,  of  appropriate  dimensions  so  that  D,  G,  Gh  are  all  compatible  with  P,  and  then 
form  the  interconnection  Pdg  =  (DPD-1  —  /3*G)Gft.  Note  that  when  forming  GGft  we  were  able 
to  employ  state  cancellations  to  obtain  a  stable  realization  with  the  same  number  of  states  as  G*,. 
It  turns  out  that  we  may  once  again  exploit  the  special  structure  of  G  and  G u  to  obtain  further 
state  cancellations  when  forming  the  interconnection  Pdg-  This  state  reduction  comes  about  by 
cancelling  uncontrollable  and/or  unobservable  modes  in  the  interconnection,  and  is  highly  desirable 
from  a  numerical  viewpoint  since  the  Hca  optimal  controller  computed  in  step  3  of  procedure  3  will 
have  the  same  number  of  states  as  our  realization  for  Pdg- 

Note  from  theorem  5  that  the  State  Space  formulae  for  the  blocks  of  Gh  and  GGft  share  the 
same  “A”  and  “B”  matrices.  In  other  words  we  have  realizations  for  each  block  of  Gh  and  GG/t 
which  take  the  form 


(G/,); 


Si] 

(nr'  t  _ 

Ai 

Si 

Ci 

D{  J 

A 

Di 

(50) 


Now  that  it  is  easy  to  show  that  since  each  block  of  G^  and  GGf,  satisfies  these  relationships,  we  may 
choose  the  block  diagonal  system  realizations  of  Gh  and  GGj,  to  satisfy  them  as  well.  Furthermore 
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we  may  scale  GGft  by  a  real  scalar  (namely  —  /?*)  and  still  preserve  these  relationships,  since  one 
may  easily  check  that  for  any  (3  G  R 


P 


A 

B  ‘ 

'  A 

B  ' 

C 

D 

[pc 

(3D 

(51) 


Thus  we  may  choose  our  realizations  for  G^  and  —  /?„GGft  to  he  of  the  form 


Gft 


Ag 

Bq 

—  ft  nn,  — 

Ag 

Bg 

Cn 

Dh\ 

.  C a 

Da 

(52) 


Now  suppose  that  our  State  Space  realization  for  DPD  1  is 


i  _ 

A  D 

bd 

cD 

Dd\ 

Then  the  interconnection  Pdg  Is  given  by 

Pdg  =  (DPD-1  —  /?*G)Gft 
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1 1  0  0\ 

where  we  have  made  use  of  the  similarity  transformation  T  =  0  I  0  J  .  Note  that  in  the  final 

\o  I  I) 

formula  for  Pdg  we  only  have  the  states  from  DPD  1  and  G/,  appearing  once  each. 

4.2.3  Fitting  D  and  G 

Finally  we  consider  the  problem  of  fitting  State  Space  system  realizations  to  frequency  response 
data.  Of  course  we  will  not  attempt  to  exactly  interpolate  the  data  at  each  frequency  point  since 
that  would  require,  in  general,  too  many  states.  Rather  we  will  attempt  to  obtain  realizations  that 
approximate  the  data  over  the  frequency  range  of  interest.  Note  that  since  we  are  concerned  with 
the  Hoo  norm  of  the  resulting  interconnection,  we  still  require  acceptable  behavior  from  the  State 
Space  fits  even  outside  the  frequency  range  we  are  especially  interested  in.  In  particular  we  require 
that  the  fits  behave  reasonably  as  the  frequency  tends  to  zero  or  infinity. 
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First  of  all  consider  the  problem  of  fitting  D(u).  This  problem  is  encountered  when  imple¬ 
menting  a  “D-K  iteration”  for  complex  p  synthesis,  and  hence  has  already  received  a  good  deal 
of  attention.  Several  methods  have  been  developed  to  tackle  this  problem,  and  although  there  is 
certainly  room  for  improvement,  they  often  work  quite  well. 

In  particular  the  //-Tools  toolbox  [BDG+91]  contains  two  algorithms  for  fitting  scalar  magnitude 
data  with  a  SISO  system.  One  algorithm,  fitmag,  uses  the  complex-cepstrum  (see  [OS75])  to 
generate  an  appropriate  phase,  and  then  uses  the  algorithm  invfreqs  from  the  signal  processing 
toolbox  [LS88],  which  uses  an  equation  error  method  to  identify  a  system  model  from  the  frequency 
response  data.  Another  //-Tools  algorithm,  fitmaglp,  uses  a  linear  programming  technique  to  find 
a  SISO  system,  which  attempts  to  fit  the  magnitude  data.  In  addition  to  these  two  algorithms  for 
fitting  magnitude  data,  there  is  the  algorithm  fitsys,  which  attempts  to  fit  given  scalar  frequency 
response  data  in  both  magnitude  and  phase  (with  a  SISO  system).  This  algorithm  is  based  heavily 
on  the  algorithm  invfreqs  (see  [LS88]  for  more  details). 

Now  one  may  easily  build  a  routine  to  fit  D{ u)  by  using  the  above  algorithms  to  fit  the  (scalar) 
elements  of  D(u>)  with  SISO  systems.  The  diagonal  elements  can  be  fitted  (in  magnitude  only) 
with  fitmag  or  fitmaglp,  and  the  non-zero  off  diagonal  elements  can  be  fitted  with  fitsys.  Putting 
these  systems  together  appropriately  we  can  form  D  £  1ZM  fitting  D(u>).  This  is  the  procedure 
we  will  use  here,  and  we  will  not  go  into  any  further  details,  since  it  is  the  same  procedure  one 
typically  uses  in  the  complex  p  synthesis  case.  Note  that  since  D  is  a  State  Space  system,  D(jw) 
tends  to  a  real  matrix  as  u  goes  to  zero  or  infinity,  and  of  course  we  have  that  M(P,K)(jw)  also 
tends  to  a  real  matrix  (for  the  same  reason).  The  following  theorem  is  easily  shown  from  results  in 
[You93] . 


Theorem  8  Given  a  matrix  M  £  R"xn  and  a  compatible 


block  structure  1C  then  we  have  that 
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where  VI,  Gl  are  defined  as 
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(55) 
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Proof:  Follows  immediately  from  the  machinery  of  theorems  5.11  and  6.1  of  [You93].  □ 


It  follows  from  this  result  that  the  optimal  scaling  D{ui),  at  zero  or  infinite  frequency,  may  be  taken 
to  be  real  as  well.  Thus  the  low/high  frequency  limits  of  D{lo)  are  amenable  to  fitting  with  a  State 
Space  system  D(ju>).  Note  also  that  all  blocks  of  D  are  square,  and  furthermore  each  block  of 
D(ju>)  may  be  taken  to  be  Hermitian  positive  definite.  Thus  we  do  not  need  to  allow  sign  changes 
in  the  (eigenvalues  of  the)  blocks  of  D(ju>)  and  hence  we  do  not  need  poles  or  zeros  on  the  imaginary 
axis.  Thus  we  may  restrict  our  fitting  routines  to  return  D  £  TZM  invertible,  with  D,D-1  £  TZCoo 
(which  is  typically  done  in  the  complex  p  synthesis  case,  for  the  same  reasons)  as  we  required  in 
section  4.2.2. 

Now  we  turn  our  attention  to  the  problem  of  fitting  jG(u>).  Once  again  we  will  tackle  this 
problem  by  fitting  the  (scalar)  elements  of  j G(u>)  with  SISO  systems,  and  using  these  to  build 
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up  G.  For  the  non-zero  off  diagonal  elements  we  can  once  again  fit  them  using  fitsys.  Note  that 
at  zero  or  infinite  frequency  theorem  8  implies  that  the  optimal  G(u)  may  be  taken  to  be  pure 
imaginary,  so  that  j G(u>)  may  be  taken  to  be  pure  real,  and  hence  amenable  to  fitting  with  a  State 
Space  system  G(jw).  For  the  diagonal  elements  of  jG(w)  we  have  additional  constraints,  which  are 
explored  in  more  detail  below. 

Note  that  we  need  to  fit  j G(u)  in  both  magnitude  and  phase.  In  particular  we  require  that  the 
phase  of  the  diagonal  elements  of  G  is  purely  imaginary  for  all  frequency.  In  order  to  see  what 
constraints  this  places  on  G  consider  first  the  conditions  under  which  a  scalar  transfer  function 
g(.s)  is  purely  real  for  all  frequency. 

Lemma  6  Consider  a  scalar  transfer  function  g  E  TZM .  Then  the  following  statements  are  equiv¬ 
alent 

1.  g(jw)  is  purely  real  for  all  u  E  K. 

2.  g~(a)  =  g(-s)  =  g(s). 

3.  g (s)  has  pole  and  zero  symmetry  about  the  imaginary  axis,  with  the  total  number  of  poles  and 
zeros  at  the  origin  being  even. 

4 .  g (s)  =  where  z(s2)  and  p(s2)  are  polynomials  in  s2. 

Proof:  We  will  prove  that  1  — »  2  — »  3  — »  4  — »  1.  Suppose  first  g(jw)  is  purely  real  for  all 
u£l.  Then 

g(j^)  =  g*(jw)  =  g(-j«) 

and  hence  g(s)  —  g(— s)  is  zero  for  all  s  =  j u>.  But  since  (g(s)  —  g(— s))  E  7 ZM  it  has  finitely  many 
zeros  and  hence  we  must  have  g(s)  —  g(— s)  =  0,  or  in  other  words  g~(s)  =  g(— s)  =  g(s).  Now 
g(s)  =  g(— s)  gives  that  whenever  s0  is  a  zero  (pole)  of  g(s),  then  — s0  is  a  zero  (pole)  of  g(s). 
Together  with  the  fact  that  for  g  E  TZM  the  poles  and  zeros  come  in  conjugate  symmetric  pairs, 
this  implies  that  we  have  pole  and  zero  symmetry  about  the  imaginary  axis.  Since  this  implies  that 
all  the  poles  and  zeros  not  at  the  origin  must  come  in  pairs  (or  fours),  it  is  easy  to  check  that  we 
must  also  have  the  total  number  of  poles  and  zeros  at  the  origin  being  even.  By  collecting  terms 
in  pairs  it  is  easy  to  check  that,  if  we  have  pole  and  zero  symmetry  about  the  imaginary  axis,  and 
the  total  number  of  poles  and  zeros  at  the  origin  is  even,  then  we  may  write  g (s)  =  Finally 

it  is  clear  that  if  g(s)  =  then  g(ju>)  is  purely  real  for  all  wER.  □ 

A  simple  extension  of  these  arguments  gives  us  the  conditions  under  which  g(s)  is  purely  imaginary 
for  all  frequency. 

Lemma  7  Consider  a  scalar  transfer  function  g  E  TZM.  Then  the  following  statements  are  equiv¬ 
alent 

1.  g(jw)  is  purely  imaginary  for  all  u  £  R. 

2.  g~(a)  =  g(-s)  =  -g(s). 

3.  g(s)  has  pole  and  zero  symmetry  about  the  imaginary  axis,  with  the  total  number  of  poles  and 
zeros  at  the  origin  being  odd. 
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4-  g(s)  =  where  z(s2)  and  p(s 2)  are  polynomials  in  s2 . 

Proof:  Again  we  prove  that  1  — »  2  — *  3  — >  4  — *■  1.  First  suppose  that  g(jw)  is  purely 
imaginary  for  all  u  G  R.  Then 

g(jw)  =  -g*(ju>)  =  -g(-jw) 

and  hence  g(s)  +  g(-s)  is  zero  for  all  s  =  j u>.  But  since  (g (s)  +  g(-s))  G  TIM  it  has  finitely  many 
zeros  and  hence  we  must  have  g(s)  +  g(-s)  =  0,  or  in  other  words  g~(s)  =  g(— s)  =  — g(s).  Now 
g (s)  —  -g(-s)  gives  that  whenever  s0  is  a  zero  (pole)  of  g(s),  then  —  s0  is  a  zero  (pole)  of  g(s). 
Together  with  the  fact  that  for  g  G  TIM  the  poles  and  zeros  come  in  conjugate  symmetric  pairs, 
this  implies  that  we  have  pole  and  zero  symmetry  about  the  imaginary  axis.  Since  this  implies  that 
all  the  poles  and  zeros  not  at  the  origin  must  come  in  pairs  (or  fours),  it  is  easy  to  check  that  we 
must  also  have  the  total  number  of  poles  and  zeros  at  the  origin  being  odd.  By  collecting  terms  in 
pairs  it  is  easy  to  check  that,  if  we  have  pole  and  zero  symmetry  about  the  imaginary  axis,  and  the 
total  number  of  poles  and  zeros  at  the  origin  is  odd,  then  we  may  write  g(s)  =  Finally  it  is 

clear  that  if  g(s)  =  s then  g(jw)  is  purely  imaginary  for  all  w£R.  □ 


Thus  we  see  that  a  purely  imaginary  transfer  function  is  obtained  by  adding  a  zero  at  the  origin 
to  a  purely  real  transfer  function.  Our  algorithm  to  fit  purely  imaginary  SISO  transfer  functions, 
fitmagreal,  which  is  used  to  fit  the  diagonal  elements  of  jG'(a/),  exploits  this  fact  and  is  described 
briefly  below. 

Suppose  that  g{ui)  is  a  diagonal  element  of  G{ u>).  Then  we  first  fit  h(u>)  =  with  a  purely 


real  SISO  transfer  function,  g.  In  order  to  do  this  we  parametrize  our  allowable  transfer  functions 
for  the  fit  as 


gO) 


■K*2) 

p(s2) 


where  the  degree  of  p{s2)  is  chosen  by  the  user,  and  the  degree  of  z(s 2)  is  chosen  so  that  g(s)  is 
strictly  proper.  Thus  we  are  left  with  the  problem  of  finding  the  coefficients  of  z  and  p  so  that 
g(jcj)  =  p\zZ^)  ^est  approximates  h(u>).  Note  that  all  the  quantities  concerned  in  this  problem 
are  real  numbers.  We  tackle  this  problem  using  an  algorithm  realfit ,  which  is  a  modification  of 
the  algorithm  fitmaglp ,  and  is  based  on  a  linear  programming  approach.  Note  that  h(u)  is  not 
constrained  in  sign,  so  that  we  wish  to  allow  sign  changes  in  g(jw).  Hence  we  do  not  place  any 
restrictions  on  the  pole/zero  locations  of  g(s),  and  in  particular  it  may  have  poles  and/or  zeros  on 
the  imaginary  axis. 

Having  obtained  g(s)  as  a  strictly  proper  transfer  function,  say 


A 
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(58) 


then  we  form  g  as 


g(s)  =  sg(s) 


[  A 

B  ' 

[CA 

CB 

(59) 


Thus  we  obtain  g  as  a  purely  imaginary  transfer  function,  with  g(jw)  approximating  jg(oj).  Note 
from  this  construction  that  the  transfer  matrix  G,  built  from  these  transfer  functions,  may  contain 
poles  and  zeros  located  on  the  imaginary  axis.  For  this  reason  the  factorization  of  G  used  in  section 
4.2.2  did  not  make  any  assumptions  about  the  pole/zero  locations  of  G. 
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Note  that  lemma  7  implies  that  any  non-zero  purely  imaginary  (proper)  transfer  function  g  £ 
1ZM  is  in  fact  strictly  proper,  and  indeed  our  construction  for  g  will  yield  a  strictly  proper  transfer 
function.  Thus  g(jw)  tends  to  zero  as  u  tends  to  infinity.  Furthermore,  assuming  g(s)  has  no  poles 
at  the  origin,  then  g(s)  has  a  zero  at  the  origin,  so  that  g(jcj)  also  tends  to  zero  as  ui  tends  to  zero. 
Now  by  earlier  arguments  we  know  that  theorem  8  implies  that  the  optimal  G(u)  may  be  taken 
to  be  purely  imaginary  at  zero  or  infinite  frequency.  But  G(u>)  is  Hermitian  and  so  the  diagonal 
elements  are  purely  real,  and  hence  the  diagonal  elements  of  the  optimal  G(u)  may  be  taken  to  be 
zero  at  zero  or  infinite  frequency.  Thus  the  high  and  low  frequency  limits  of  g(jw)  are  compatible 
with  the  data  it  is  attempting  to  fit  (i.  e. ,  both  j£f(w)  and  g(jw)  tend  to  zero  at  high/low  frequency). 

Finally  we  note  that  since  the  routine  fitsys  can  be  used  to  fit  a  SISO  transfer  function  in  both 
magnitude  and  phase,  we  could  apply  this  algorithm  directly  to  the  problem  of  fitting  the  diagonal 
elements  of  G{ui),  rather  than  the  approach  taken  above.  Note  however  that  the  above  approach 
enforces  the  phase  constraints  for  all  frequency,  and  enforces  the  correct  behavior  at  the  high  and 
low  frequency  limits.  We  would  not  obtain  these  guarantees  by  using  fitsys.  Furthermore  by  simply 
applying  fitsys  we  would  not  be  exploiting  all  the  information  that  we  have  about  the  structure 
of  g  (as  is  done  above),  and  in  fact  we  found  that  the  fitsys  routine  performed  poorly  for  many 
problems  when  fitmagreal  was  able  to  obtain  quite  good  fits. 

4.2.4  Some  Practical  Considerations 

These  procedures  have  been  implemented  in  software  as  Matlab  functions  (m-files)  “rmufit”  and 
“rmuflp.”  Together  with  the  mixed  //,  analysis  software,  “rmu”  (see  [YND92]),  they  allow  the  user 
to  implement  the  “ D,G-K  iteration”  described  in  the  preceding  subsections,  and  hence  to  design 
mixed  //  synthesis  controllers. 

It  should  be  noted  that  for  complex  /t  problems  then  the  “D.G-K  iteration”  reduces  exactly 
to  the  familiar  “ D-K  iteration”  for  complex  p  synthesis.  Thus  many  of  the  techniques  required 
to  obtain  a  practical  design  from  a  complex  //  synthesis  procedure  carry  over  into  the  mixed  case 
as  well.  There  are  many  issues  to  be  considered,  and  we  will  not  go  into  any  of  the  details  here, 
except  to  make  the  following  comments.  The  procedures  developed  in  this  paper  did  not  make  any 
restrictions  on  the  block  structure,  1C,  of  the  uncertainties  in  the  underlying  robustness  problem.  In 
particular  repeated  scalar  blocks  (real  or  complex)  were  allowed.  It  should  be  noted  however  that 
for  real  or  complex  fi  problems  the  presence  of  repeated  scalar  uncertainty  blocks  gives  rise  to  full 
blocks  in  the  scaling  matrices  (D(oj),  G(u?)).  Since  the  present  approach  fits  element  by  element, 
then  this  may  well  require  a  lot  of  of  states,  and  this  approach  may  become  impractical. 

Fortunately  for  many  problems  we  do  not  need  repeated  scalar  uncertainty  blocks,  and  so  we 
only  need  diagonal  scaling  matrices.  In  fact  the  currently  available  software  in  the  //-Tools  toolbox 
(for  complex  n  synthesis)  only  handles  this  case.  For  this  reason  most  of  the  effort  in  our  approach 
has  been  to  develop  a  method  which  will  work  well  for  that  case  (hence  the  emphasis  on  developing 
good  routines  for  fitting  the  diagonal  elements  of  D(u>)  and  jG(w)).  A  practical  solution  to  the 
more  general  case  is  a  subject  of  ongoing  research. 

5  Examples  and  Applications 

Finally  we  conclude  the  paper  with  a  look  at  two  simple  design  problems  involving  real  parametric 
uncertainty.  In  both  cases  it  will  be  seen  that  a  simple  application  of  complex  //,  synthesis  techniques 
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is  not  adequate,  but  the  new  mixed  fi  synthesis  procedure  performs  well. 

Optimal  Gain  Margin  Problem 

The  problem  of  designing  a  controller  to  maximize  the  gain  margin  of  a  SISO  plant  can  be  solved 
analytically  (see  section  11.3  of  [DFT92]).  This  question  can  (almost)  be  restated  as  a  special  case 
of  a  mixed  fi  synthesis  problem,  so  that  it  provides  us  with  a  simple  benchmark  test  for  the  “ D,G-K 
iteration”  procedure. 

Consider  the  feedback  interconnection  in  figure  4.  The  perturbed  plant  P  is  given  by 

p{s)  =  (r=rs) (1 + n 

where  Sr  E  K  is  a  real  scalar  uncertainty.  This  is  a  loopshaping  set-up,  and  fits  into  the  standard 
mixed  /r  synthesis  framework  of  figure  1.  If  we  choose  the  performance  weights  for  sensitivity  and 
complementary  sensitivity  as 


WSl  =  Ws 3  =  WTl  =  WT2  =  0.01 

then  the  performance  requirements  for  this  problem  are  negligible.  Thus  we  find  that  this  robust 
performance  problem  is  effectively  a  robust  stability  problem,  and  is  in  fact  (approximately)  the 
optimal  gain  margin  problem. 


Figure  4:  Interconnection  for  mixed  n  synthesis  examples 

Note  that  the  nominal  plant  is  non-minimum  phase  and  unstable.  It  can  be  stabilized  simply 
with  unity  negative  feedback,  and  the  Nyquist  plot  is  shown  in  figure  5  (for  Sr  =  0).  In  addition  the 
mixed  and  complex  \i  plots  across  frequency  for  the  closed  loop  system  are  also  shown.  Note  that 
in  this  example  the  mixed  n  upper  and  lower  bounds  were  tight,  so  that  only  the  upper  bounds  are 
shown.  For  a  plant  of  the  form 


1  —  as 
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with  a  >  1,  then  the  optimal  gain  margin  controller  is  given  by  [Doy92] 


K  =  - 


1  +  fis 
s  +  f3 


where  /3  j  -.  The  results  for  the  optimal  gain  margin  controller  are  shown  in  figure  6,  and  it  can 
be  seen  that  this  controller  flattens  out  the  Nyquist  curve  so  as  to  substantially  improve  the  gain 
margin.  As  a  result  the  peak  value  across  frequency  of  mixed  fi  is  substantially  improved  (from  6.0 
to  3.35),  at  the  expense  of  complex  fi. 


Nyquist  Plot  Mu  Plots  For  Unity  Feedback 


1 )  Complex  Mu  2)  Mixed  Mu 


Figure  5:  Nyquist  and  fi  plots  for  unity  (negative)  feedback  controller  (gain  margin  problem) 
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Figure  6:  Nyquist  and  fi  plots  for  optimal  gain  margin  controller  (gain  margin  problem) 

The  standard  “ D-K  iteration”  procedure  for  complex  fi  synthesis  was  applied  to  this  problem, 
and  the  results  are  shown  in  figure  7.  Note  that  this  controller  gives  a  slight  improvement  over 
unity  feedback  for  both  mixed  and  complex  fi  (from  6.0  to  5.57).  However  the  Nyquist  plot  is  not 
greatly  modified  (other  than  the  roll  off  at  high  frequency),  so  that  the  gain  margin  is  not  greatly 
improved.  Applying  the  mixed  fi  synthesis  procedure  to  this  problem  results  in  the  plots  in  figure 
8.  Note  that  this  controller  achieves  about  the  same  level  for  mixed  fi  as  the  optimal  (gain  margin) 
controller.  Once  again  this  is  achieved  by  flattening  out  the  Nyquist  curve  so  as  to  improve  the  gain 
margin  at  the  expense  of  complex  fi.  For  this  simple  example  we  see  that  the  mixed  fi  synthesis 
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procedure  does  exactly  the  required  thing,  and  shapes  the  frequency  response  so  as  to  minimize 
the  peak  value  across  frequency  of  mixed  /x. 


Nyquist  Plot  Mu  Plots  For  Complex  Mu  Controller 


Figure  7:  Nyquist  and  /x  plots  for  complex  /x  synthesis  controller  (gain  margin  problem) 


ISTyqiiist  Plot  Mu  Plots  For  Mixed  Mu  Controller 


Figure  8:  Nyquist  and  /x  plots  for  mixed  /x  synthesis  controller  (gain  margin  problem) 


A  Modal  Example 

The  previous  example  involved  only  one  uncertain  real  parameter,  so  that  it  was  possible  to  exploit 
our  knowledge  of  the  problem,  and  find  a  good  controller  design  by  another  method  for  the  purposes 
of  comparison.  We  found  that  the  mixed  /x  synthesis  routine  was  able  to  find  a  comparable  design. 
Of  course  the  main  point  in  developing  the  mixed  /x  synthesis  machinery  was  to  have  a  systematic 
design  procedure,  capable  of  dealing  with  more  complicated  problems,  where  alternative  approaches 
are  no  longer  available.  We  conclude  the  paper  with  a  synthesis  example  for  a  problem  with  two 
uncertain  real  parameters. 

Consider  once  again  the  interconnection  in  figure  4.  As  in  the  gain  margin  problem  we  make 
this  essentially  a  robust  stability  problem  by  choosing  the  weights  as 

WSl  =  Ws2  =  WTl  =  WT2  =  0.01 
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so  that  the  performance  requirements  are  negligible.  Now  consider  an  uncertain  plant  of  the  form 

p,  \  -2C^  +  o;2(l  +  ^) 

W  s*  +  2(us  +  ^(l  +  8r2) 

where  8[,82  6  R  are  uncertain  real  parameters.  This  is  a  second-order  transfer  function  intended 
to  represent  a  lightly  damped  mode  in  a  flexible  structure.  This  can  be  put  in  standard  //  synthesis 
form  by  rewriting  as  P  as  in  figure  9  where  A  =  diag(tf[,  S2)  and 
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For  this  example  we  choose  u>  =  1  and  (  =  0.2,  so  that  the  mode  is  lightly  damped,  with  a  pair  of 


Figure  9:  Modal  example  uncertain  plant 

complex  conjugate  poles  at  —0.2  ±  0.9798J  and  a  pair  of  complex  conjugate  zeros  at  0.2  ±  0.9798j 
(for  =  S2  =  0). 

Consider  for  a  moment  what  happens  to  these  pole/zero  locations  as  we  start  to  move  8[,82 
away  from  zero.  If  we  restrict  8[,82  G  R  then  we  get  the  motion  indicated  by  the  arrows  in  figure 
10,  so  that  initially  the  poles  stay  in  the  left  half-plane  and  the  zeros  stay  in  the  right  half-plane. 
However  if  we  allow  8[,82  G  C  then  we  get  motion  inside  the  indicated  circles.  In  this  case  we  find 
that  the  poles  and  zeros  can  cross  the  imaginary  axis  for  much  smaller  values  of  the  uncertainty. 
Thus  we  would  expect  to  see  a  fairly  large  gap  between  mixed  and  complex  fi  for  this  type  of 
problem. 

A  complex  /i  synthesis  controller  was  designed  for  this  problem  and  the  complex  and  mixed  //, 
plots  are  shown  in  the  right-hand  plot  of  figure  11.  Once  again  only  the  upper  bounds  are  shown. 
This  controller  achieves  a  peak  value  for  complex  /i  of  about  2.6,  with  mixed  //  being  about  the 
same.  The  results  for  the  mixed  /r  synthesis  controller  are  shown  in  the  left-hand  plot  of  figure  11. 
This  controller  does  not  greatly  affect  the  peak  value  of  complex  yz,  but  reduces  the  peak  value  of 
mixed  \i  to  about  1.1.  Thus  we  do  indeed  find  that  there  is  a  large  gap  between  mixed  and  complex 
l-i  for  this  problem,  and  the  mixed  //,  synthesis  procedure  is  able  to  exploit  that  fact  to  significantly 
improve  on  the  complex  /z  synthesis  controller. 
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Figure  10:  Pole/zero  motion  for  the  modal  example 


Complex  TVIvx  Synthesis  Controller 


IVIix-ed  Mu  Synthesis  Controller 


Figure  11:  Complex  and  mixed  n  plots  for  the  modal  example 


Of  course  the  two  designs  presented  here  are  only  for  simple  example  problems,  where  the 
controller  is  chosen  to  optimize  the  robust  stability  margin  with  respect  to  some  uncertain  real 
parameters.  A  practical  control  design  for  a  real  problem  involves  a  great  many  other  issues.  In 
particular  one  would  need  to  select  weights  so  as  to  reflect  a  more  suitable  performance  criterion, 
rather  than  merely  optimizing  the  stability  margin  at  the  expense  of  everything  else  (see  [YA]). 
Nevertheless  these  examples  do  illustrate  the  workings  of  the  new  method,  and  the  initial  results 
are  encouraging.  A  further  example  may  be  found  in  [YA],  where  this  design  technique  is  compared 
to  classical  methods  for  a  robust  performance  problem  with  real  parameter  uncertainty. 

6  Conclusion 

The  general  mixed  \x  synthesis  problem  has  been  considered,  and  a  “D,G-I(  iteration”  procedure 
developed  to  tackle  the  problem.  This  procedure  relies  heavily  on  the  theoretical  properties  of 
fi  and  its  bounds,  and  has  guaranteed  convergence  to  a  local  minimum  of  the  problem  (which  is 
nonconvex).  This  yields  a  stabilizing  controller  which  attempts  to  minimize  the  peak  value  across 
frequency  of  mixed  n,  and  hence  provide  the  best  possible  robust  performance  in  the  face  of  mixed 


31 


(real  and  complex)  uncertainties.  The  controller  synthesis  procedure  is  applicable  to  general  mixed 
g  robust  performance  problems,  although  further  work  is  required  to  develop  a  practical  solution 
to  the  repeated  parameter  case.  This  procedure  has  been  implemented  in  software  and  applied  to  a 
number  of  example  problems  involving  real  parametric  uncertainty,  where  it  was  able  to  synthesize 
a  controller  yielding  a  lower  value  for  mixed  g  (and  hence  better  robust  performance)  than  could 
be  obtained  with  standard  complex  g  synthesis  techniques. 
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